Signatures of collective flow in high multiplicity pp collisions. 
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A blast- wave parametrization, including a full set of hadronic resonances, is used to model a small 
system, with total particle multiplicity comparable to the one measured in the high-multiplicity pp 
collisions at the LHC. Calculations are preformed for three cases: with negligible, regular and strong 
radial flow on the blast-wave hypersurface. We investigate the effects of flow on inclusive pr spectra 
as well as on ID and 3D femtoscopic radii for pions. Special emphasis is put on the role of pions from 
resonance decays, fn particular we show that they magnify the flow effects present in the blast-wave 
stage and significantly influence the shape of the correlation functions. A specific observable, the 
Rout/Rfide ratio is proposed as a sensitive probe of the collective effects. Model results for the high 
multiplicity pp collisions, for scenarios with small and large radial flow are compared. 
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I. INTRODUCTION 

The pp program at the LHC has successfully started, 
and has already produced a large set of minimum-bias 
measurements from the experiments [IHl]. Of particu- 
lar interest is the fact that the multiplicity range mea- 
sured at y/s = 7 TeV overlaps with the one measured in 
peripheral CuCu collisions at ultra-relativistic energies 
of ^Wn = 62 and 200 GeV at RHIC. This allows for 
comparison of observables between p+p and heavy-ion 
collisions at the same measured final state multiplicity 
dNch/dq. In heavy- ion collisions, dN c h/drj is the scal- 
ing variable for many soft physics observables, includ- 
ing the inclusive particle pt spectra, pion femtoscopic 
radii [9(, strange particle yields etc. It is therefore inter- 
esting whether a high-multiplicity pp collision resembles 
a heavy-ion one. 

This question has become even more relevant with 
the recent observation of the long-range r\ ridge in high 
multiplicity pp collisions by the CMS experiment [lOj. 
There are attempts to explain the phenomenon in the 
Color Glass Condensate framework [Tljj or as an elliptic 
flow [1^ . Both of these critically depend on the obser- 
vation of radial flow (or space-momentum correlation') in 
such collisions. A recent study with EPOS model [13[ 
suggests that such flow has actually been seen in the 
ALICE pion femtoscopy data Similar topics have 
been investigated in the framework of the rescattering 
model [14]. A comparison of pion femtoscopic radii in pp 
and heavy-ion collisions at lower energy y/s = 200 GeV 
have been carried out by the STAR experiment [l5| . Scal- 
ing between radii in pp and heavy-ion collisions has been 
observed vs. the pair momentum, suggesting that there 
might be a common physics mechanism, i.e. radial flow, 
driving them in both cases. 

The hydrodynamic scenario is well applicable to heavy- 
ion collisions at RHIC. It produces specific space- 
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momentum correlation patterns, which are commonly re- 
ferred to as flow. The system created in the heavy-ion 
collision expands rapidly outwards, showing a strong ra- 
dial flow, which is observed in the modification of the 
single particle inclusive pt spectra shape. The spectrum 
is an observable depending only on the momenta of the 
particles, so its connection to space-time can only be 
indirectly inferred. To directly access space-time infor- 
mation we employ femtoscopic techniques. The fall of 
the "femtoscopic radii" (sometimes called "HBT radii") 
with particle's m<r can be interpreted as the decrease of 
"lengths of homogeneity" , a direct consequence of radial 
and longitudinal flow [la]. The simultaneous description 
of the pt spectra and the femtoscopic radii dependence 
on rriT in heavy-ion collisions can be achieved with the 
"blast-wave" parametrization In this work we 

test what results such an approach brings when applied 
to small systems. We characterize the size of the system 
by the average charged particle multiplicity per unit of 
pscudorapidity that it produces. This means that sim- 
ple "blast-wave" models cannot be used, as abundances 
of each particle species are free parameters there. In 
contrast the THERMINATDR model [5(| has both necessary 
features: it can use the traditional "blast- wave" hyper- 
surface and provides absolute particle multiplicities. 

The usage of the THERMINATOR makes it possible to 
test the impact of the strongly decaying resonances. As 
it was shown in [l9j , their introduction increases the ap- 
parent size of the system produced in heavy-ion collisions 
by approximately 1 fm and leads to the development of 
long-range tails in the emission function of pions. These 
effects were important but not dominant in the AuAu 
system (with the average size of 6 fm), but they may 
dominate in a pp collision where we expect the initial 
system to have a size of the order of 1 fm. 

The paper is organized as follows. In Section |H] we 
briefly introduce the THERMINATDR model and its blast- 
wave variant. In Section IIIII we describe the choice 
of parameters for the three simulations that were per- 
formed. In Section IIVI we describe the results for the 
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inclusive transverse momentum pt spectra and the 1- 
dimensional femtoscopic simulations and explain the fea- 
tures observed there. In Section|V]we show the results of 
the 3D femtoscopic analysis and propose the most robust 
and sensitive observables that can be used to search for 
collectivity in the pp collisions at large multiplicities, e.g. 
at the LHC. In Section l*VTl we discuss the origin of the 
effects seen in previous Sections. 

II. THERMINATOR MODEL AND FEMTOSCOPIC 
FORMALISM 

In this work we use the THERMINATOR model (20|, more 
specifically its variant that uses the blast- wave freeze-out 



hypersurface with the possibility to adjust the space-time 
correlation at the freeze-out hypersurface. For a detailed 
description of the model, its full theoretical basis and 
formalism see [19]. In short, average per-event parti- 
cle abundances are calculated from the chemical model, 
where the temperature T and bariochemical potential \ib 
are model parameters, and isospin fii and strange ps po- 
tential are calculated from conservation laws. For each 
event the number of particles of each type (we use 381 
particle types taken from PDG (21]) is randomly gener- 
ated from a Poissonian distribution. Then each particle 
is put on the freeze-out hypersurface according to the 
generalized blast-wave prescription [19J : 
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where a, t, /9 max , are parameters of the model. p max is 
the transverse size of the system, r is the proper time 
of emission, and a determines the transverse position (p) 
vs. time (t) correlation at the hypersurface. Of particular 
importance is the form of the radial flow velocity, which 
depends on particle emission position: 

Pi Pmax / n 

v ±- = ; — / » (2) 

where vt is a parameter which controls the amount of 
flow in the system. Large values of vt give small flow, 
small values of vt give large flow, growing semi-linearly 
from at the center of the source to the maximum value 
of 1/(1 + i>t) at the edge, y, pt and ip are the components 
of the momentum of the particle (rapidity, transverse mo- 
mentum and azimuthal angle respectively), while an, p 
and (f> are the corresponding coordinates of the emission 
point. As the last step, all unstable particles are prop- 
agated and decayed, in cascades if necessary, until only 
stable particles remain. The emission point of the par- 
ticle is either its creation point on the hypersurface (we 
call such particles "primordial") or a decay point of its 
parent resonance. 

Final state particle rescatterings after the emission are 
not included. We have estimated that for the particle 
densities in this work, each particle would rescatter on 
the average 0.8 times. Most of the rescattering points 
would be very close to the original emission point, so the 
inclusion of rescattering would not change the emission 
shapes significantly. 

Further in the work we show results of the analy- 
sis of the identical pion femtoscopic correlations (some- 



times called "Bose-Einstein correlations" or "HBT cor- 
relations"). All the analyzes were performed using the 
standard "two-particle" method: pions from the same 
event were combined into pairs, and assigned a weight 
equal to their wave-function squared, to form the sig- 
nal. The pair wave function contained only the relevant 
symmctrization; final state coulomb and strong interac- 
tions were ignored as for pions they are small and are 
not important for the determination of the femtoscopic 
radii. Then, pions from different events were combined 
into pairs to form the background. The correlation func- 
tion was constructed by dividing the signal by the back- 
ground. It was then fitted with various functional forms, 
which will be described in detail in Section IIVI For an 
extensive description of the two-particle method and the 
formalism see 1191 . 



III. SIMULATION DESCRIPTION 

We performed the simulations with the blast-wave 
model with the following parameters: transverse size of 
the freeze-out hypersurface p m ax = 1.5 fm, proper time 
of freeze-out To = 2.0 fm/c, freeze-out temperature Tf = 
165.6 MeV, bariochemical potential ps = 28.5 MeV. 
The freeze-out hypersurface was slightly sloped in the 
p — t plane (a = —0.2). The three simulations differed 
only in the value of the vt parameter; the "no-flow" had 
vt = 10, the "regular flow": vt — 1 and the "strong 
flow": vt = 0.5. The resulting flow velocity profiles can 
be seen in Fig. [T] as dashed lines. The common radial flow 
present in the model results in strong space-momentum 
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FIG. 1. (Color on-line) Average common transverse flow of 
primordial pions (see Eq.®) as a function of the distance to 
the central axis of the source p. "No flow" results are shown 
as blue circles, "regular flow" as red triangles, "strong flow" 
as green squares. The input model flow velocity profiles (see 
Eq. p])) are shown as dashed lines. 
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FIG. 2. (Color on-line) Single particle inclusive mr spectra 
for pions (green squares), kaons (red triangles) and protons 
(blue circles). Upper panels a) and b) show simulations with 
"strong flow", lower panels c) and d) with "no flow". Left 
panels a) and c) show all particles, right panels b) and d) 
only primordial particles. 



p — pt correlations of the produced particles, which is 
reflected in the the common radial flow: 



P 



(3) 



where averaging is done over particles. The results are 
shown in Fig. Q] as symbols. The plot illustrates that 
the "no flow" scenario has very little p — pt correlations, 
while in the "strong flow" scenario this correlation is sig- 
nificant. For each of the scenarios we have simulated 2M 
events. 



IV. RESULTS 

Fixing the model parameters, as described above, fixes 
the average absolute multiplicity of particles per event. 



TABLE I. Slope parameters of the single particle spectra 
shown in Fig. in MeV. 
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FIG. 3. (Color on-line) Kaon to pion (panel a)) and proton 
to pion (panel b)) ratios vs. pt- Green dashed lines show 
primordial particles for "no flow", red long-dashed show all 
particles for "no flow", blue dot-dashed show primordial for 
"strong flow", violet full line show all for "strong flow". 



The resulting (dN c h /drj) for the "strong flow" sample was 
8.3. This is about twice the average multiplicity reported 
by ALICE collaboration at </s = 900 GeV and is above 
the averag e multiplicity for ^fs = 7 TeV 

The impact that the radial flow has on particle spec- 
tra is shown in Fig. [2] Only for the simplest case of "no 
flow" and primordial particles only (shown in panel d)) 
the spectra for all particle types are exponential. We fit 
the distributions with an exponential and extract their 
slope, shown in Tab. |U Even in the "no flow" case the 
simple addition of strongly decaying resonances obscures 
the picture, as seen in panel c). Resonances tend to pop- 
ulate the lower tut regions more abundantly than higher 
rriT- There are relatively more resonances decaying into 
pions and protons than into kaons. As a result, the slope 
parameter is higher for pions and protons and lower for 
kaons. For primordial particles in "strong flow" in panel 
b) we see, as expected, that the thermal "exponential" 
shapes are modified by radial flow: pions get a "concave" 
shape, while kaons and protons develop a positive curva- 
ture. All are less steep than in the "no flow" case, re- 
sulting in the larger slope parameter. Adding resonance 
decay products, shown in panel a), decreases the slope 
parameter by 50 MeV, so any attempt to extract the 
flow velocity and temperature from the slope parameters 
must fully take into account the resonance contribution. 
In Fig.[3]the kaon to pion and proton to pion ratios vs. pt 
are shown. For K/tt both flow and addition of resonances 
decrease the ratio. This comes from the larger contribu- 
tion of resonance decays to pions. In contrast, for p/ir 
only flow decreases the ratio, as the relative contributions 
of resonances to both pions and protons are similar. In 



4 



summary, the particle spectra retain the general charac- 
teristics associated with radial flow even when combined 
with resonance decays in a small system, but resonance 
decays do alter the quantitative estimates of tempera- 
ture significantly. In addition our simple model does not 
include the hard processes which will contribute to all 
particle's spectra at larger and complicate the picture 
further. Therefore making strong conclusions about ra- 
dial flow from particle spectra alone is not trivial and is 
model dependent. 

Another observable which was extensively used in con- 
junction with pt spectra as a signature of collective flow 
is the pion femtoscopic radius of the system, more specif- 
ically its monotonic decrease with growing pair momen- 
tum kx = \Pt,i +Pt,2\ /2. The first measurement of 
such a dependence at the LHC was reported in Q. The 
minimum-bias measurement shows no significant depen- 
dence, however the combined multiplicity vs. kx analysis 
shows that the behavior of "low" and "high" multiplic- 
ity events is different. The "low" multiplicity (M) events 
show an increase of pion femtoscopic radius Rinv at low 
kx and then a fall at high kx- The "high" M events on 
the other hand suggest a decrease of the Ri nv radius with 
kx, although it is also consistent with a flat dependence. 
New results from the LHC, especially from collisions at 
the higher energy of y/s = 7 TeV will deliver results with 
better precision and will show if the trend continues at 
higher multiplicities. 

We have obtained the pion femtoscopic correlation 
functions using the two-particle method described in 
Section [TTJ We divide the sample into six kx bins: 
(0.12,0.2), (0.2,0.3), (0.3,0.4), (0.4,0.5), (0.5,0.6) and 
(0.6,0.7) GeV/c. We begin by analyzing the 1- 
dimensional correlation functions. Following the experi- 
ments [3,0]; we n t them with two functional forms, the 
Gaussian: 



fa(qinv) = 1 + Xcxp(-Rf nv 2 qf n 



,), 



and the exponential: 



fE(qinv) = 1 + \eiq>(~Rf nv q inv ) 



(4) 



(5) 



An example of the fits is shown in Fig. SJ We see a 
behavior similar to the one reported by the LHC experi- 
ments: the correlation function is better described by an 
exponential. The Gaussian fit is worse, but it character- 
izes the overall width of the correlation function. 

In Fig. [5] we plot the results of the 1-dimensional anal- 
ysis vs. kx for the three scenarios of flow. The Gaussian 
and properly scaled exponential fit results are consistent. 
In the "no flow" scenario the Ri nv is growing with kx, 
similar to low multiplicity collisions at the LHC. On the 
other hand, as flow develops, the Rf nv becomes flat vs. 
kx for "regular flow" and has a negative slope for "strong 
flow" . In other words, from the comparison of the Rf nv 
vs. kx trends between data and the model, we observe a 
development of radial flow with the increase of multiplic- 
ity in the pp collisions at the LHC. It would be interesting 
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FIG. 4. (Color on-line) 1-dimensional pion femtoscopic 
correlation function for the "strong flow" case, kr bin of 
(0.2,0.3) GeV/c. Triangles are the results of the simulation, 
dashed line is a Gaussian fit, dotted line is the exponential 
fit. 
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FIG. 5. (Color on-line) Overall 1-dimensional femtoscopic 
radius Ri nv of the system. "No flow" are shown as blue circles, 
"regular flow" as red triangles, "strong flow" as green squares 
(the latter two shifted slightly for clarity). Gaussian radii 
Rf nv are shown as closed symbols, exponential radii Rf nv as 
open symbols. 



to see if the data at higher multiplicities at y/s = 7 TeV 
continue the trend and develop an even stronger kx de- 
pendence. The results are intriguing enough to attempt 
a more in-depth explanation of their origin. We aim to 
propose more strict tests, which would make the conclu- 
sions from such comparisons less model dependent. 



V. 3D BERTSCH-PRATT ANALYSIS 

We have performed the 3D analysis of the simu- 
lated correlation functions. We have used the standard 
Bertsch-Pratt decomposition of the relative momentum 
into the "long" (along the beam axis), the "out" (along 
the pair transverse momentum), and the "side" (per- 
pendicular to the other two) directions. We have used 
the Longitudinally Co-Moving System (LCMS) reference 
frame, where the longitudinal momentum of the pair van- 
ishes. We have analyzed two classes of pions: the "pri- 
mordial" only, and "all" particles. The latter includes 



5 



a - 

O 1 



■a. 



o regular, k y (0.3, 0.4) GeV/c 
|q |<30 MeV 

'^other' 



0.1 0.2 0.3 0.4 0.5 

q out (GeV/c) 




0.1 0.2 0.3 0.4 0.5 

q.„ (GeV/c) 

side 




q (GeV/c) 

"org v ; 

FIG. 6. (Color on-line) Projections of the 3-dimensional 
correlation function for the "strong flow" simulation for 
kr{0.3, 0.4) GeV/c. Dashed lines show fits to 1-d projections. 
Panel a) shows the out projection with an exponential fit. 
Panel b) shows side with Gaussian, panel c) shows long with 
exponential. To plot the projections, other components were 
integrated in the range |g ther| < 0.03 GeV/c. 



the "primordial" plus those from strongly decaying reso- 
nances. In Fig. |5]we show the correlation function, sepa- 
rately in out, side and long directions. We attempted to 
fit both the Id projections and the full 3d function with 
combinations of Gaussian, exponential and Lorentzian 
shapes. These three were chosen because they have an 
analytic relation between correlation function and pair 
emission function forms. We found that (a) for primor- 
dial particles a 3D Gaussian is the best fit, (b) for all 
particles the out and long direction is best fitted with an 
exponential, while side remains Gaussian. This is true 
both for Id projections (shown in the figure) as well as 
for a full 3d fit. Even though the fit in Fig.|5]is best out of 
three forms tried, it is still not perfect. Nevertheless, the 
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FIG. 7. (Color on-line) kr dependence of the 3-dimensional 
radii for the primordial particles. Blue circles are "no flow" 
simulations, red triangles are "regular flow", green squares 
are "strong flow". Panel a) shows R^ ut , b) shows Rfide, c) 
shows Rf ong , d) shows Rout/R?ide ratio. 



width of the function is adequately estimated, which is 
enough for the purpose of this paper. The investigation 
of the details of the shape of the correlation functions we 
leave for future studies. 

Based on the study of correlation function shapes, for 
primordial particles a Gaussian 3d ellipsoid, with differ- 
ent sizes in all directions is a reasonable description of 
the correlation: 



R 



long 



2 \ 
llonq) 



out j R-side an d R-fong are Gaussian radii in out, 



where 

side and long directions respectively, A is the strength 
of the correlation. For all particles we have used the 
following fit functional form: 

Cs(g) = ^ + ^ ex P(~Rout\loiLt\-Rf l de <£ide-RFong\Uong\), 

(7) 

where Rg Ut and Rf ong are exponential radii in out and 
side directions. Similarly to the 1-dimensional case, the 
radii obtained from Gaussian and exponential fits are not 
directly comparable, as they are parameters of different 
functional forms. 

In Fig [7] we show the Gaussian femtoscopic radii ob- 
tained from the analysis of "primordial" particles. We 
concentrate on the transverse dynamics. For the R^ ut 



6 



o 



LU o 

DC 



o 



DC 



_ 

£ 3 



LU .O 

DC 



DC 



DC 1 



no flow 
regular 
strong 



0.2 



0.4 

k T (GeV/c) 



0.6 



FIG. 8. (Color on-line) kr dependence of the 3-dimensional 
radii for all particles. Blue circles are "no flow" simulations, 
red triangles are "regular flow", green squares are "strong 
flow". Panel a) show Rout, b) shows Rfide, c) shows Rf ong , 
d) shows Ro Ut /Rf ide ratio. 



radius we observe an expected behavior: the "no flow" 
scenario shows a flat kx dependence, the slope grows with 
increasing flow and is significant for "strong flow" . In 
addition, the size gets smaller with increasing flow; the 
well-knows mechanism of the "lengths of homogeneity" 
is at work here. Rf ide dependence is weakly changed by 
flow; stronger flow creates a small slope. The change in 
behavior is best visible in the Rf ut /Rf ide ratio. For "no 
flow" it is large, at least 1.0-1.5, for "strong flow" it is 



falling, and smaller, below 1.0. The Rf u tl Rfide ratio has 
an advantage that it is, to a limited degree, independent 
of scale. The absolute values of the transverse radii are 
directly depending on the assumed system size (the p m ax 
parameter), while the ratio is less dependent on it, and 
its qualitative behavior should not change as we change 
Pmax in a range which would produce multiplicities ob- 
served in pp collisions. 

For primordial particles we have identified the 
R'out/Rfide ra tio as a sensitive, scale-independent, probe 
of the amount of collectivity in the system. However, 
the majority of pions are coming from resonance decays, 
which strongly influence spatial and temporal character- 
istics of the emission process. It is therefore crucial to 
investigate how they modify the signal present for "pri- 
mordial" particles. In Fig. [3] we show the fit results for 
all pions. We stress that inclusion of resonances changes 
the emission patterns visibly, which is reflected in the 
change of the fit functional form. R„ ut for "no flow" is 
larger than for "primordial" and develops a kr slope. In 
contrast, the "strong flow" R„ ut grows little: R® ut grows 
less as flow grows. That is surprising, as the "length of 
homogeneity" argument seems to hold for particles from 
resonances, even though they do not come directly from 
a "collective" medium. Also the negative slope of the 
Rout dependence on kx alone is not a good signature of 
collectivity as such behavior also develops for "no flow" . 
Rf ide shows the behavior opposite to Rf ut . For "no flow" 
it is similar to "primordial" case, while for "strong flow" 
it grows. As a consequence, again in striking contrast 
to the behavior of R„ ut , Rf ide grows more as radial flow 
increases. Again the Rout /Rfide ratio is the most sen- 
sitive to these changes, and for all particles it shows an 
even more dramatic difference between flow scenarios: 
the value for the ratio is up to a factor of 1.5 larger for 
"no flow" as compared to "strong flow" . In summary, 
we show that a small value of Routl Rfide ratio, close to 
unity, is a signal of emission of all particles from a collec- 
tive medium. The value is monotonically dependent on 
the strength of flow: the stronger the flow, the lower the 
Rout/ Rfide ratio. In addition, the inclusion of resonance 
decay products magnifies this flow effect. 



VI. ORIGINS OF FEMTOSCOPIC RADII 
BEHAVIOR 

The identification of the Routl Rfide ratio as a sensi- 
tive probe of collective particle emission in systems re- 
sembling pp collisions at the LHC has consequences for 
experiments. It allows to test a hypothesis of a creation 
of a collective system in such collisions, which produces 
large particle multiplicities. In addition, the femtoscopic 
pion analysis is practically a "minimum-bias" measure- 
ment and does not require advanced triggering or accep- 
tance corrections. It can therefore be performed quickly, 



with the data samples that the LHC experiments already 
have on tape. It is therefore critical to understand why 
such effect arises and how general their emergence is; in 
other words how model dependent is the link between 
Rout/ Rfde ratio and system collectivity. 

The separation between two particles in the out direc- 
tion is: 

rout = xi - x 2 + /J(*2 - h) (8) 

where x = ppt cos(</>— ip)/pT are the transverse emission 
point coordinates of the first and the second particle, 
projected on the pair velocity, t\ and t% are the emission 
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FIG. 9. (Color on-line) Correlation between space and time in particle emission from the "no flow" model. Panels a), b) and 
c) show "outwards" direction (x from Eq. ((5|): projection on particles' velocity. Panels d), e) and f) show "sidewards" (y from 
Eq. ([9])) direction: projection perpendicular to the velocity. Panels a) and d) show primordial particles, b) and e) resonance 
products, c) and f) all particles. The graphs overlaid on the plots show the mean and spread of the distribution in a time bin. 
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FIG. 10. (Color on-line) Correlation between space and time in particle emission from the "strong flow" model. Panels a), b) 
and c) show "outwards" direction: projection on particles' velocity. Panels d) , e) and f ) show "sidewards" direction: projection 
perpendicular to the velocity. Panels a) and d) show primordial particles, b) and e) resonance products, c) and f) all particles. 
The graphs overlaid on the plots show the mean and spread of the distribution in a time bin. 



times of the particles and /3 is the pair velocity (we ne- 



glect the differences between particles velocities since in 



8 



order to be femtoscopically correlated their relative ve- 
locity must be small). We order the particles in a pair 
in such a way that t% > t±, the third term in the for- 
mula reflects the fact that the first particle has already 
traveled some distance before the second was born. The 
definition of the side direction is that the pair velocity in 
this direction vanishes, so the formula for the separation 
is simpler: 

Tside =2/1-2/2 (9) 

where y = (ppt sm((j)—ip))/pT are the transverse emission 
point coordinates of the first and second particle, now 
perpendicular to the pair velocity, x coordinates are on 
the x axis of the top panels of Figs. and [TU1 y are on 
the x axis of the bottom panels and t is on the y axis of 
all panels. The R„ u f and R%de femtoscopic radii have 
the meaning of variances of out and side distributions 
respectively (conventionally divided by v2)- They are 
proportional to a combination of the single particle space- 
time distribution moments: 

Ron? « C» + fat - C(X, fit) (10) 
Rfide ~ Oi, (11) 

where a x , a y and at are the variances of the x, y and t 
distributions respectively and C(x, fit) is the term that 
accounts for a possible space-time correlation in the emis- 
sion process. We note that C contributes to the size with 
a negative sign: it means that positive x — t correlation 
decreases the out radius [HI, [23j| . 

In Fig. |9] the space-time characteristics of the emission 
process for the "no flow" scenario is shown. On panels a) 
and d) are primordial particles. In this case both out and 
side directions show similar size, which does not depend 
on time. Also, the mean emission points are not shifted. 
Panels b) and e) show the behavior for resonance decay 
products. There is a small correlation between mean x 
and time. Also, out and side spreads grow with time, but 
only slightly. As a consequence, for all particles shown in 
panels c) and f), the out emission point is weakly corre- 
lated with time and the side spread grows only slightly. 

In Fig. [10] we present a similar study for the "strong 
flow" scenario. We concentrate on the differences be- 
tween this scenario and the "no flow" one. For primor- 
dial particles, the out mean position is shifted by flow. 
The spread is slightly smaller than side; hence the smaller 
than one ratio for green points in panel a) of Fig. [7] For 
resonance products we see more striking differences. The 
out mean position is strongly positively correlated with 
time. The spread also grows with time, but less than for 
the side, where the growth is more pronounced than for 
the "no flow" case. As a result, for all particles, the out 
mean position is positively correlated with time, while 
the side spread grows stronger than in the "no flow" case. 

We can now understand the dramatic change of the 
Rout/Rfide r& tio between "no flow" and "strong flow" 
scenarios. For "no flow" , a x and a y are similar, as can 
be seen in panels c) and f ) of Fig. [§] But R^ ut is further 



increased by at and the x — t correlation is too small to 
counter-balance it. As a result Rout/Rfide 1S l ar S e - I n 
short, resonance decays introduce significant spread in 
emission times, which acts to increase Rf ut - The fact 
that resonance decay points are distributed according to 
the exponential decay time multiplied by velocity also 
explains why the shape of the correlation in the out di- 
rection changes dramatically, while it stays Gaussian in 
side, where the velocity vanishes. 

In the "strong flow" scenario we note that even for pri- 
mordial particles, seen in panel a) and d) of Fig. [TQl a x 
is smaller than a y - the effect of flow of primordial pions. 
The effect is further increased by resonances products. 
Because of the existence of common radial flow in the 
system, when resonances are born on the freeze-out hy- 
persurface, they are already flying to the outside of the 
system. Since they decay after some time, the emission 
points of the daughter pions will, on average, be further 
apart compared to the "no flow" case. This is nicely il- 
lustrated as the strong growth of the a y with time on 
panel e) of Fig. [TQl The result is a significant increase 
of R% de , as compared to the "no flow" case. One ex- 
pects the same argument to hold for a x , and to some 
extent it is seen in panel b) of Fig. [TO] However, the 
growth is not as strong as for a y - this is because the 
resonance emission points are also governed by flow, so 
their parent a x is smaller than a yi just like for primordial 
pions on panels a) and d). But more importantly, radial 
flow of resonances also results in much stronger (as com- 
pared to "no flow" case) positive x — t correlation, seen 
most clearly in panel b). This correlation also dominates 
when we consider both primordial and resonance daugh- 
ter pions in panel c). This will further decrease the R„ ut 
parameter, due to the C(x, fit) term in Eq. [JOJ As a re- 
sult, the R^ut I Rfide ratio is small. In short, radial flow of 
primordial pions and resonances produces several effects 
which all combine together to bring the R® u tl R%d& ratio 
down, even below unity if the flow is strong enough. The 
stronger the flow, the lower the ratio. 

We also now understand why 1-dimensional Ri nv radii 
dependence on kr changes as flow develops. Rinv is 
approximately equal to the quadratic average of 7 Rout, 
Rside and i?/ onff , where the 7 factor in front of R ou t comes 
from the boost (with pair velocity) from the LCMS to 
Pair Rest Frame, where Rinv is determined. For "no 
flow" case R ou t is flat with kr in LCMS, so the growth 
of the boost as kx increases drives Ri nv up. For cases 
with flow, the growth of the boost is countered by the 
fall of Rout resulting from flow. And since the other two 
radii also decrease with kx, the slope of the Ri nv depen- 
dence on kr becomes negative. 

As a final note we stress that the model employed here 
is a simple hydrodynamics-inspired parametrization. We 
have shown that one of its particular features - an p — px 
correlation, also known as radial flow, produces specific 
signatures in single particle spectra as well as ID and 
3D pion femtoscopic radii. However it does not address 
the question of where such correlations come from. One 
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possibility is that they come from a "medium" that is de- 
scribed by hydrodynamic equations, as suggested in [l3j . 
but there may be other mechanisms that produce such 
correlations. 



VII. CONCLUSIONS 

We have performed simulations with the THERMINATOR 
model, employing emission of primordial particles from a 
blast-wave hypersurface, and including the propagation 
and decay of a full set of hadronic resonances. The pa- 
rameters of the calculations have been adjusted, so that 
they produce multiplicities observed in high multiplicity 
pp collisions at LHC energies. Three sets of parameters 
have been used, which differ only by the strength of the 
radial flow (the p — pt correlation) in the system. 

The radial flow modifies the particle spectra in the ex- 
pected way, creating a concave shape for pions and posi- 
tive curvature for kaons and protons. However, addition 
of resonance decays significantly modifies these shapes, 
influencing the extracted slope parameter and making a 
direct extraction of flow velocity from spectra measure- 
ments alone complicated and model dependent. 

Results of 1-dimensional femtoscopic analysis show 
that with no flow, the ID Rf nv pion femtoscopic radius 
is expected to grow with pair momentum kr- In con- 
trast, as flow is increased, the dependence changes to a 



decrease with kx- Therefore, within the frame of this 
model the results of ALICE collaboration on multiplic- 
ity vs. kx dependence can be interpreted as signature of 
the development of radial flow with increasing per-event 
particle multiplicity. 

3D pion femtoscopic radii show that the shape of the 
correlation is exponential in the out and long direc- 
tions due to contributions from resonance decay prod- 
ucts, while in side the shape remains Gaussian. Radial 
flow (or lack of it) for primordial pions and resonances has 
opposite effect on Rf ut and R% de radii. The first grows 
less as flow develops, which is a result of a combination 
of two effects: (a) decrease of lengths of homogeneity 
for primordial pions and resonances and (b) strong x — t 
correlations coming from resonances' flow-ordered veloc- 
ities. Rj ide grows more as flow develops, the effect of the 
preferential propagation of the resonances to the outside 
of the source in the presence of flow. We propose to use 
the Ro U t/Rfide ratio as a probe of the amount of radial 
flow (e.g. velocities ordering) in the small, resonance- 
dominated system, produced in the high multiplicity col- 
lisions at the LHC. Any analysis aiming to search for 
collective effects in such systems should combine at least 
the three measurements: identified particle pt spectra, 
ID and 3D femtoscopic radii. 
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